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CN . Abstract 

pi • Mapping the central 5 pc region of the nearby radio galaxy NGC 

, 6251 with a 0.2 pc resolution using VLBI at two radio frequencies, 

5 GHz and 15 GHz, we have found the sub-pc-scale counter jet 

(N ■ 

for the first time in this radio galaxy. This discovery allows us to 
investigate the jet acceleration based on the relativistic beaming 
i model (Ghisellini et al. 1993). 

> ■ 
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. 1 Introduction 

. The genesis of powerful radio jets from active galactic nuclei is one of 

I the long standing problems in astrophysics (Bridle & Perley 1984). Al- 

' though global morphological properties give us very important informa- 

■ tion, very inner regions in the radio jets also provide hints to understand- 

ing the genesis of radio jets. In order to investigate radio galaxies at the 
O ' sub-milli arcsecond angular resolution, we have performed new high- 

-4— > ' resolution VLBI observations of NGC 6251 using HALCA (Hirabayashi 

S : et al. 1998). We use a distance to NGC 6251, 94.4 Mpc (for a Hubble 

^ I constant Hq = 75 km s~^ Mpc~^). Note that 1 mas (milli arcsecond) 

corresponds to 0.48 pc at this distance. 

2 Observations 



NGC 6251 was observed at 5 GHz using VSOP on 30 April 1998 and 
at 15 GHz using VLBA on 2 June 1998. Details of the observations 
are summarized in Table 1. In order to perform beam-size-matched 
comparison between 5 GHz and 15 GHz, we restored the two images 
with a same spatial resolution of 0.50 x 0.50 mas. 
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Table 1: Observations. 



Date 


V 

(GHz) 


Stations 


Peak Int. 
(Jy/beam) 


RMS noise 
(m Jy/beam) 


DR' 


1998 Apr. 30 


5 


VLBA, HALCA, EB'' 


0.13 


0.50 


260 


1998 Jun. 02 


15 


VLBA 


0.34 


0.25 


1400 



^DR : Dynamic range; ^EB : Effelsberg 100-m telescope 



Figure 1: The images of NGC 6251 at 5 GHz (a) and at 15 GHz (b). 
The images are rotated clockwise on the sky by 28°. The spectral index 
map is also shown (c). Here the spectral index a is defined as Si, (x u"'. 



3 Results 

Our final maps at 5 GHz and 15 GHz are shown in Figure la and lb, 

respectively. Although the secondary peak is seen at both frequencies, 
its angular distance from the 15 GHz brightest peak is larger by 0.3 mas 
(0.14 pc) than that at 5 GHz. This difference may be attributed to the 
different opacity toward the central engine at 5 GHz and 15 GHz. 

Here we apply two new methods for the registration using correla- 
tion functions between intensity profiles and between wiggling patterns 
of the radio jet at the two frequencies. The correlation function of inten- 
sity profiles exhibits the maximum correlation coefficient 2± 0.98 at an 
offset of 0.30 ± 0.11 mas, indicating that the peak at 5 GHz is shifted 
to the jet direction from that at 15 GHz. Analysis of wiggling patterns 
gives an offset of 0.60 it 0.36 mas with a correlation coefficient of ~0.66. 
Since these two offsets were obtained using the two independent meth- 
ods, we adopt a weighted average of the two offsets. Thus the best offset 
of 0.42 ± 0.14 mas is obtained. Using this offset, we register the two 
images and obtain a spectral index image of the radio jet as shown in 
Figure Ic. It can be seen optically thin component in the opposite side 
of the jet. We conclude that this is the real counter jet. 
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Figure 2: A schematic illustration of the spherical- absorption model (a), 
and spectral index variation along the radio jet of NGC 6251 (b). 



4 Discussions 

We assume for simplicity that the core is surrounded by a plasma sphere 
with a radius of a and the radio emission from the inner part of the jets 
(i.e., both the jet and the counter jet) suffer from the free- free absorp- 
tion. The approaching jet escapes from the plasma sphere at a projected 
distance of x = a x sin^jet and thus the spectral index becomes to be 
intrinsic here. The counter jet suffers the effect of free-free absorp- 
tion until X = —a. It is also noted that the path length is longest at 
X = —a X sin^jet. As shown in Figure 2a, we define the following pro- 
jected distances; Xj^t = ^peak = a x sin^jet and X^et = o- Then we 
are able to estimate O-^^t = sin~^ (x^)' adopt the 15 GHz bright- 
est peak as the core, because the optical depth at 15 GHz is enough 
to be small in the case of free-free absorption. In Figure 2b, we show 
the observed spectral index variation along the radio jet. We estimate 
-'^peak ~ 0.24 mas, Xjet ~ 0.41 mas, and Xcjet ~ 0.80 mas. Then we 
obtain O^^t — 31°. 

Adopting the so-called Doppler beaming model, the jet to the counter 
jet intensity ratio R can be written as a function of projected distance x. 



= [ i-g::(gc::fc(g j ' ' > where /3jet(x) = vUx)/c. In Figure 
3a, we measured i? as a function of x. It is shown that R is estimated 



-a{x) 
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Figure 3: The jet/counter jet flux density ratio R as a function of the 

projected distance from the brightest peak at 15 GHz for the 5GHz 
data {filled triangle) and for the 15 GHz data {filled circles) (a), and jet 
velocity /3jet as a function of projected distance from the brightest peak 
at 15 GHz (b). 



to be systematically larger at 5 GHz than those at 15 GHz. This is 
probably due to stronger absorption at 5 GHz than at 15 GHz. It is 
shown that R at 15 GHz increases from 1.7 at 0.85 mas (0.4 pc) to 10 at 
1.25 mas (0.6 pc) with projected distance from the core. If this increase 
is just caused by the Doppler beaming, it is suggested that /9jet increase 
with distance, because it cannot be seen that ^jet varies significantly. 
Figure 3b shows that the jet is accelerated from /?jet ~ 0.1 at 0.4 pc 
to /3jct ~ 0.5 at 0.6 pc. This provides the first direct evidence for the 
acceleration at the sub-pc-scale radio jet. 
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